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Abstract. This paper reexamines univariate reduction from a toric geometric 
point of view. We begin by constructing a binomial variant of the w-resultant and 
then retailor the generalized characteristic polynomial to fully exploit sparsity in 
the monomial structure of any given polynomial system. We thus obtain a fast new 
algorithm for univariate reduction and a better understanding of the underlying 
projections. As a corollary, we show that a refinement of Hilbert's Tenth Problem 
is decidable within single-exponential time. We also show how certain multisym- 
metric functions of the roots of polynomial systems can be calculated with sparse 
resultants. 



1. Introduction 

We give a new approach for combatting degeneracy problems which occur when 
reducing large polynomial systems to univariate polynomials. Taking a positive atti- 
tude, we will actually make use of these degeneracies to better understand polynomial 
system solving. We do this by applying new observations involving toric varieties to 
establish faster, more reliable algorithms for univariate reduction. Our techniques 
provide an intrinsic geometric setting for such reductions and fully exploit the spar- 
sity of any polynomial system specified by its monomial term structure. 

Algebraically reducing problems involving polynomial systems to simpler questions 
involving a single univariate polynomial is a technique which has been known for over 
a century, and has been applied with increasing efficiency in computational algebra 
over the last two decades. (For example, see ||Laz81| , jKE^ , |FEI)3| , |GHMF95|| and the 



references therein.) However, there are difficulties with this method which still have 
not been completely addressed. Three such problems are the following: 

A: Since univariate reduction corresponds to projecting the solution space onto a 
projective line, the fiber over a point with finite coordinate might contain roots 
at infinity. 
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B: The univariate polynomial one reduces to might be identically zero, due to (1) 
degeneracies in the underlying projection, or (2) the presence of infinitely many 
roots. 

C: Univariate reduction is still too slow for many problems of interest. 
We circumvent A by the technique of toric laminations. By constructing a gener- 
alization of Canny's generalized characteristic polynomial (GCP) ||Can90| , |Shu93| — 
the toric (or sparse) GCP — we can dispose of B2. Our toric geometric approach also 
helps us better understand Bi. These two new techniques are contained in theorems 
1-3 below. 

Although problem B2 was solved earlier in ||Ren87| , |Can90|| , the toric GCP greatly 
improves the complexity bounds given there. In particular, we address problem C by 
using the sparse resultant [|GKZ90|, Pj^C95|| throughout our development. Fast methods 



Emi94 



for polynomial system solving were derived via the sparse u-resultant in 
li)C95|| but problems A and B were not addressed there. Our approach to univariate 
reduction thus unifies the GCP with the the sparse it-resultant. Furthermore 
generalize both techniques to certain toric varieties |[KKMS73| , |l^'ul93| , [Stu95 



we 



One advantage of a toric variety setting is reducing the problem of exact root 
counting for polynomial systems, almost completely, to convex geometry. Generically 
sharp upper bounds on the number of (complex) isolated roots, in terms of mixed 
volume, first appeared just over two decades ago ||Kus75|| . These bounds have since 
been generalized to algebraically closed fields | |Dan78 , ^ioj97a[ | and to various subsets 
of affine space (other than the algebraic torus) ||Kho78| ~ HS96 , |Roj97a|] . However, the 



following algorithm appears to be the first example of a convex geometric approach 
to counting the exact number of roots. (We will assume henceforth that all our 
polynomials and roots are considered over an algebraically closed ground field K.) 

Theorem 1. Let F{xi, . . . ,x„) be an nxn polynomial system with support E := 
{El, . . . ,En) such that Ai{E) > 0. Also let be the sum of the convex hulls of 
the Ei, and pick a G Z'^\{0} such that the segment [0,a] is not parallel to any 
facet of Pe- Define Tia{u^,uJ) := Res£;^(F, it+ + tt_a;") where {u±} is a pair of 
algebraically independent indeterminates and Ea := {E,{0,a}). Finally, let e± be 
the lowest exponent of u± occurring in any monomial of iTa- Then F has exactly 
Af := A4{E) —£+ — £_ isolated roots (counting multiplicities) in (K*)'^ , provided 
J\f<oo. 

In the above, Res*(-) and A4{-) respectively denote the sparse resultant and mixed 
volume |[EC95| , PGH96| , [Roj97a|| . The key contributions of the above theorem are (a) 
fibering the solution space by maps more general than coordinate projections, and 
(b) making explicit use of roots at toric infinity when doing univariate reduction. 
We may also calculate the number of distinct roots in a similar way. 

Corollary 1. Following the notation and assumptions of theorem 1, let 7r^(M_(_, it_) 
be the square- free part of7Ta{u-^,u^). Also let e'^ be the lowest exponent ofu± occuring 
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in TT^. Then F has exactly A/'': = deg(7r^) — — e'_ roots in (K*)" , provided Af' <oo 
and the projection of roots {(e{K*)^ \ F{Q=0} — > {C"} is injective. 

This corollary is proved alongside theorem 1 in section ^ 

By instead calculating a particular coefficient of tTq, we can actually compute cer- 
tain multisymmetric functions of the roots of F and thus do more than just count the 
number of isolated roots. Admittedly, such a computation can be rather difficult, but 
recent advances in interpolation techniques, e.g., ||Zip93|] , are making this approach 



increasingly feasible. So let us at least give explicit formulae for the coefficients of tTq 
in terms of multisymmetric functions. 

Definition 1. For any gM"\{0}, let Ef he the set of points of Ei having minimal 
inner product with w . Also, let E^ := {Ef , . . . , i?^) and, when w eQ"'\{0} , letp^jE 
IT' he the first lattice point not equal to the origin encountered along the ray generated 
hy w. Finally, let L{E,a), the ambiguity locus, he the suhvariety of Tpj^\{K*)"' 
corresponding to the union of all faces of having an inner normal perpendicular 
to the segment [0,a]. 

Remark 1. The toric variety (over K), Tp, corresponding to a polytope P is detailed 
KSZ92 , b'ul93 , LjKZ94 , Stu95 , Prloj97a| ]. Polynomial roots in a toric compactifica- 



m 



tion are descrihed (in a manner closest to our present framework) in | )Roj97a 



Corollary 2. Following the notation of theorem 1, assume instead that a is any point 
in Z"\{0} and <oo. Then is a homogeneous polynomial of degree A4{E) 

and, for any c/G {0, . . . ,Ai{E) — — £_}, the coefficient of nl^^^u^J^^^^"'^"'^'^ in Tia 
is precisely 

EH COT, 

where C G K* is independent of d, the sum ranges over all cardinality d multisets 
- C Zp, and Zp is the multiset of roots of F in (K*)" U L{E,a). 

Furthermore, if e+ = e- = Q then C is precisely 



J]Ress.(F)- 



■Pw-a 



where the product ranges over all inner facet directions of Pp having negative inner 
product with a. ■ 

Multisets are simply sets where repeated elements are allowed, and whose subsets 
have repetitions appropriately restricted ||GKP94|| . For example, if Zp contains — 1 



as a root of multiplicity 5, then no subset of Zp can have more than five — I's. 
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So in essence, the coefficients of vTa (when suitably normahzed) are just elementary 
symmetric functions of a-monomials of certain roots of F. 

The above corollary follows easily from theorem 1 and the Pedersen-Sturmfels 
product formula |pS93|| .p] Moreover, since the sparse resultant is invariant under 
muliplication of the individual polynomials by a monomial, we can compare the 
coefficients of Ha and 7r_a to obtain the following beautiful identity.^ 

Corollary 3. Following the notation and assumptions of corollary^, £+ = e_ = =^ 
where the product ranges over all inner facet directions of Pe- B 



Recall that a ridge is a polytope face of codimension 2 | |Zie95 |. Although the "bad" 



case of theorem 1 (A/'=oo) is an annoyance, it is a small annoyance and can actually 
be made use of in certain cases. 

Theorem 2. Following the notation and assumptions of theorem 1, M =oo ■^=^ F 
has a root lying in L{E,a) or infinitely many roots in {K*Y . Also, within the space 
of all systems with support contained in E, the F with M = oo form a subvariety of 
codimension >2. 

Thus, if the roots of F avoid the ambiguity locus, F has infinitely many roots in 
(K*)" <^=^ A/' = oo. In particular, it follows directly from the definition of L{E,a) 
that codimL{E , a) = 2, under the assumptions of theorem |^. (Since L{E, a) then 
corresponds to a finite set of ridges [[t:'ul93| , [Roj97a|| .) 



We can completely avoid amiguity loci by replacing the binomial m+ + in our 
preceding results with another specially chosen indeterminate polynomial. This is 
detailed further in ||Roj97c|| and, to a certain extent, in the result we quote below. 



However it is worth noting that the roots of F always avoid L(E, a) when n<3 and 
E is the support of F (cf. lemma |ip. 

Remark 2. The existence of infinitely many roots of F at toric infinity and the 
existence of a root of F within the ambiguity locus are independent in general (cf. 
section Thus it is worthwhile to know something about L{E, a). 

Remark 3. A correspondence useful for visualizing theorems 1 and 2 is to identify 
toric infinity (Tp\{K*)"^ ) with the boundary of the polytope P (cf. definition Com- 
puting TTa then amounts to splitting OPe into two signed halves, much like cracking 



^We also need the fact that the Pedersen-Sturmfels formula, o riginally stated only over C, remains 



true over a general algebraically closed field. This is proved in | Roj97d |. 

^ Alicia Dickenstein pointed out the if = C of the above formula (due to Cattani, Dickenstein 
and Sturmfcls) to the author at the conference where this paper was presented. Shortly after, the 
author realized that this formula could also follow easily from corollary ||. 
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an egg-shell. In particular, the half with sign ± consists precisely of those faces of 
Pe having an inner normal with ± inner product with a. The ambiguity locus is then 
precisely the common boundary of these two halves. We will see later that Ha can be 
viewed as a coordinatization of a particular rational map Tp^^ F]^. Thus, provided 
no roots lie in the ambiguity locus, e± is the sum of the intersection numbers of the 
roots of F which lie in the ± half. 

A general way to deal with the existence of infinitely many roots of F within a 
given toric compactification is the following result announced below. 



Theorem 3. ||Roj97c|| Following the notation of theorem 1, suppose g{x) '■=YleeA''^eX^ 



where the Ue are algebraically independent indeterminates and AdU'' is nonempty 
and finite. Suppose further that V := (Pi, . . . , P„) is an n-tuple of integral polytopes 
in R" containing the support of F , and let D be an irreducible fill of V. Setting 
^* iJ2eeDi^'^ M^fl--'^]); define H{s;u) :=Res(^E,A){F — sF* , g) , where s is a new 
indeterminate. Finally, considering H as a polynomial in s with coefficients in K[u], 
let J-'a{u) be the coefficient of the lowest power term ofTi. Then Ta is a homogeneous 
polynomial with the following properties: 

1. J-'a is divisible by 5'(C); for any root (K*)"' of F. 

2. If V is compatible with Conv(y4), then J-'a has degree A4(V) and to every ir- 
reducible factor of Ta there naturally corresponds a root of F in Tp, where 

We callli. a toric (or sparse) generalized characteristic polynomial for {F,A). ■ 

The combinatorial portion of the above theorem can be explained as follows: We 
say that D fills V iS D := {Di, . . . , Dn) satisfies Di C P^ for all i G [l..n] and 
M{D) = M{V) |[Roj94| , |RW96|| . An irreducible fill is then simply a fill which is 
minimal with respect to n-tuple containment. Finding such a D amounts to a combi- 
natorial preprocessing step which need only be done once for a given set of problems, 
provided E is fixed. Compatibility is defined in section |^ and the toric GCP is briefly 
compared to the original GCP in section ^ Combined with the toric variety version 
of Bernshtein's theorem ||Roj97a[|, we see that the above theorem gives us a resultant- 



based method to work with the isolated roots of a sparse polynomial system, even in 
the presence of infinitely many roots. Furthermore, the roots defined by Ta which do 
not correspond to any isolated root of F can be usefully interpreted as contributions 
from the excess components. 

We apply our main theorems to some simple examples in section ^ Our results 
are also very useful for root counting over M, Q, and even Z. For example, one can 
combine standard univariate techniques such as Sturm sequences ||GLRR94| , [Roy95 



with our preceding results to quickly count the number of real roots ||Roj97b|| . Better 



still, we can go to an even smaller ring: In section ^ we will give a short proof 
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(under some mild hypotheses) of the following refinement of Hilbert's Tenth Problem 
||Mat93|| . 



Theorem 4. ||Roi97d|| Consider, for any nonnegative integer d, the following collec- 



tion of Diophantine systems: 

Hilb(d): Multivariate polynomial systems with integer coefficients and d-dimensional 
complex solution set. 

Then there is an algorithm which, given any instance of Hilb(O), finds all integer 
solutions, or certifies that there are none, within single- exponential time. ■ 

All of the above results stem from an observation on the vanishing of the sparse 
resultant. In particular, one basic property of the sparse resultant is that 

(*) {F, fn+i) has a root in (iT*)" =^ Res^(F, /„+i) =0 

(provided E = {E, En+i) and fn+i is a polynomial over K with support contained in 
En+i C Zi"). However, the converse assertion is not always true. (The polynomial 
system (x + j/+l,a; + y + 2,x + ?/ + 3) with E = {0, ei, 62}^ gives a simple counter- 
example reducing to the vanishing of a 3 x 3 determinant.) The correct statement in 
which both implications hold seems to be known only folklorically, while the unmixed 
case Ei = - ■ ■ = En+i = A (over C) is contained in various recent works, e.g., ||KSZ92|, 



GKZ94|| . So we will make use of the following more general result. 



Vanishing Theorem for Resultants. [ [lloj97d| ] Suppose fi is a polynomial over K 



with support contained in i?jCZ" for all iG [l..n + 1]. Then, provided 

M{E^..... ,En+i)>0 

for some i G [l.-'n.], 

n+1 

Res^(/i,... ,/„+i)=0^ f|Pp^(/i,Conv(i?,))^0, 

1=1 

n+1 

where E:={Ei, . . . ,-E„+i) and Pe-= Y1 Conv(£'j). ■ 

i=l 

In the above, Vp{fi, Pi) is the toric effective divisor of Tp corresponding to (/«, Pi) (cf. 
definition This result provides a geometric analogue, over a general algebraically 
closed field, of the product formula for the sparse resultant [[PS93|] . 
We close this introduction with an important note: 

In addition to the applications we have presented, our approach to univariate reduc- 
tion provides a first step toward an intrinsic method for root finding in the algebraic 
torus (K*)"'. Although one could in principle reduce to various coordinate subspaces 
to perform "sparse" elimination theory in K"^, such an approach runs into compli- 
cations with intersection multiplicities on the coordinate hyperplanes and is really 
not intrinsic to K"'. We propose a more intrinsic (and efficient) approach to affine 
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elimination theory by defining a new resultant operator — the affine sparse resultant 
— in 



[Roj97e|| . Following this route, all of our main theorems generalize easily to 
root finding in affine space minus an arbitrary union of coordinate hyperplanes. This 
is pursued in greater depth in ||Roj974 |Roj97e|| and the complexity bound of theorem 



H is further refined in the latter work. 

2. Related Work 

We point out a few important recent approaches to elimination theory which can 
profit from our results. 

From an applied angle, our observations on degeneracies and handling polynomial 
systems with infinitely many roots nicely complement the work of Emiris and Canny 



EC95|] . In particular, their sparse resultant based algorithms for polynomial system 
solving can now be made to work even when problem B occurs. Also, an added bene- 
fit of working torically (as opposed to the classical approach of working in projective 
space) is the increased efficiency of the sparse resultant: the resulting matrix calcu- 
lations (for polynomial system solving) are much smaller and faster. In particular. 



whereas it was remarked in ||Can90|| that Grobner basis methods are likely to be faster 



than the GCP for sparse polynomial systems, the toric GCP appears to be far more 
competitive in such a comparison. 

From a more theoretical point of view, our focus on univariate reduction is a useful 
addition to Sturmfels' foundational work on sparse elimination theory in the large 
Stu93|| . (We also point out that [[Stu93|| provides a wonderfully clear introduction to 



some of the toric variety techniques we refer to.) 

Going further into the intersection of algebraic geometry and complexity theory, 
there is a beautiful new approach to elimination theory founded by the school of 
Heintz, et. al. | |PK95| , |GHMP9^ . Our results fill in some toric geometry missing from 



their more algebraic framework. In particular, we propose that toric laminations, 
i.e., "fibering" Tp^ by toric projections (instead of fibering (X*)" by coordinate 
projections), can enhance their methods. For example, although their formulation in 
terms of straight-line programs is extremely general, it appears that their complexity 
bounds can be significantly improved with our techniques in the case of polynomial 
systems which are "sparse" in the specified supports sense we use here. This is 



explored further in ||Roj97d |. 



3. Examples 

We will give a 2 x 2 example of our first two theorems. 

3.1. Root Counting with Sparse Resultants. Let n = 2 and consider the bivari- 
ate polynomial system F\ = (x^ + — 1, x"^ + i/^ — 1). Also let Ei and E2 respectively 
be the supports {O, (3, 0), (0, 4)} and {O, (4, 0), (0, 5)}. Clearly then the Newton 
polygons are two triangles and the segment [O, (1, 1)] is not parallel to any edge of 
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the quadrilateral Pe '■= Conv(0, (7,0), (0,9), (4,4)). So we can set a := (1, 1) and 
apply theorem 1 to count the roots of F. 

We are now ready to compute the sparse resultant ResEa{F,u+ + u^xy), where 
Ea: = {El, £'2,{0,(1,1)}) and the coefficients of F have been replaced by algebraically 
independent indeterminates.0 To do this, we can use Maple to calculate the Sylvester 
(2 X 1) resultants of the pairs {f2,u+ + U-xy) and (/i,ii+ + u^xy), eliminating the 
variable y. Let us respectively call these resultants /_i and /_2. Using the Sylvester 
resultant one last time on (/_i,/_2) to eliminate x, we arrive at a multiple of the 
sparse resultant we are looking for. (This fact follows easily from the characterization 
of the sparse resultant as the (nonzero) polynomial of lowest total degree satisfying 
property (*) from the introduction [ ^tu93| , [Roj97cl|| .) 



The preceding computations, modulo typing, takes only fractions of a second on 
a Sun 4m Computeserver. From | [Stu94| ] we know that the sparse resultant itself 
should have degree M{Ei,E2) + M{Ei,{0,a}) + M{E2,{0,a}) = 16 + 7 + 9 = 
32. (This calculation is easily done by hand via the general formula A4{Ei, E2) = 
Area(Conv(£;i + E2)) - Area(Conv(^i)) - Area(Conv(£;2)) jScEMI , PGH96|| .) So by 



factoring with Maple (which takes just under 2 seconds) we can isolate the sparse 
resultant and specialize coefficients to obtain 

7r(i_i)(M+,n_) = ResE^x^ + y^ -l,x^ + y^ - l,u+ + U-xy) 



So £+ = 7, £_ = 0, and theorem 1 tells us that F has exactly 16 — 7 = 9 roots, 
counting multiplicities, in (C*)^. Furthermore, a one-line Maple calculation shows 
that 7r(i 1) is irreducible over Q. So corollary 1 (and basic Galois theory) tells us that 
the multiplicity of every root of F in (C*)^ is 1. 

The roots of F can also be counted via Grobner basis methods but the toric 
geometry in this example is rather nice to observe. In particular, the coordinate 
cross {xy = 0} is naturally embedded in and corresponds precisely to the two 
lower-left edges of Pe (cf. remark ^ and definition |^). Furthermore, the ambiguity 
locus consists of two distinct points which, due to the fact that is not identically 
(cf. theorem 2), the roots of F have luckily avoided. Also £+ = 7 =^ F has precisely 
7 roots, counting multiplicities, on the coordinate cross. The latter fact can also be 
easily checked with a little commutative algebra. 

We have included timing information solely for illustrative purposes. In particular, 
our calculations can be sped up tremendously with suitably specialized code. 



■'The published version of this paper contains an erroneous computation (due to a software 
error) of this resultant. This was observed by John Dalbec at the conference where this paper was 
presented, and the calculation below is due to him. The author is deeply indebted to Professor 
Dalbec's astute observation. 
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4. A Refinement of Hilbert's Tenth Problem 



We will present a concise proof of theorem 4 under the following three hypotheses: 

1. The instances of Hilb(O) considered are n x n polynomial systems. 

2. The complex roots of an instance must have all coordinates nonzero. 

3. No instance can have a complex root at toric infinity. 

Note that (1) and (2) are of a combinatorial nature: they restrict the number and 
dimensions of the supports, as well as their intersections with the coordinate sub- 
spaces. On the other hand, (3) is of a more algebraic nature and, when the supports 
are fixed, is easily shown to fail only on a codimension 1 subvariety of polynomial 
systems, e.g., ||Roj974 Main Theorem 2]. 

The above hypotheses are present mainly for technical reasons and are removed in 
| Roj97d |. They may also be removed with other recent techniques |PPR94 , PK95 , 



|GHMP95|1 , but through our toric techniques it is possible to obtain a general com 
plexity bound with near quadratic dependence on the mixed volume ||Roj97c| 



This 



would be the best asymptotic bound to date for solving Diophantine systems. 
Proof of Theorem 4: First we outline our algorithm: 



Step 1: Following the notation of theorem 3, let A be the vertices of a standard n-simplex 

in and V the n-tuple of Newton polytopes of F. Form the polynomial JF4. 
Step 2: Using the Lenstra-Lenstra-Lovasz algorithm ||LLL82|| , factor J-'a over Z. 
Step 3: Any integral root {zi, . . . ,Zn) G of F corresponds precisely to an integral 
factor of J^A of the form uq + ZiUg^ + ■ ■ ■ + z^Ue^. 
Step 1 runs in time polynomial in the mixed volume, applying hypotheses (1) and 
(3), and the results of ||Stu93| , [Emi96|| . In fact, we are merely computing a "sparse" 
variant of the m- resultant in Step 1. To prove this (and the fact that this sparse 
M-resultant is not identically 0) one uses slightly more general versions of theorems 1 



and 2 |[Roj97d|| . One also needs the fact that hypotheses (1) and (3) together imply 
that the complex zero set of F in (K*)^ is zero-dimensional, but this is an immediate 
consequence of [ [Roj97a] , Corollary 3]. We also point out that working more generally 
in terms of J-'a aids in later removing hypothesis (3) ||Roj97d[ . 

To see that Step 2 runs within time single-exponential in the input size, one need 
only observe that (a) the coefficient growth in the sparse w-resultant calculation can be 
reasonably bounded, and (b) the LLL-algorithm runs within time single-exponential 
in n and polynomial in the bit-sizes of the coefficients of F. Part (a) is detailed in 
Ro,i97d|| , and part (b) follows easily from the univariate case covered in |[LLL82 . 



That we can terminate as described in Step 3 follows immediately from hypothesis 
(2) and the Vanishing Theorem for Resultants. ■ 

We conclude this section with some brief notes on related work. 

First, we point out the rather surprising fact that restructuring Hilbert's Tenth 
Problem as the union [J^^g Hilb((i) seems to be new. (In particular, Hilbert's Tenth 
Problem was originally stated as deciding the existence of a single integral root for 
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one polynomial in several variables.) Although the results of ||Jon82|| imply that 
Hilb(8) is undecidable, nothing seems to be known about Hilb((i) for < < 8. 
We find this shocking, considering the vast effort within arithmetic geometry to use 
complex geometric invariants in Diophantine problems, e.g., Falting's theorem and 
the deep conjectures of Lang and Vojta ||CS86|| . 

We also point out that Hilb(O) being a single-exponential problem may not be so 
new: Teresa Krick and Luis-Miguel Pardo-Vasallo have pointed out to the author 
that theorem 4 can also be readily derived from [[PKQSI , pHMP95[] . (Indeed, such an 
argument could have also been derived directly from |[LLL82 | and [ [Ren87|| in 1987.) 
Pardo-Vasallo has also pointed out that Marie-Frangoise Roy lectured in 1995 on a 



similar result, presumably derived from 



Nevertheless, optimal bounds are 



open territory; not to mention the 6- long gap in the status of Hilb((i). 

5. ToRic Varieties and Sparse Resultants 



Our notation is a slight variation of that used in |Pul93|| , and is described at 
greater length in ||Roj97a]| . We will assume the reader to be familiar with normal 
fans of polytopes and the construction of a toric variety from a fan | Ful93 , G^KZ94 . 
we will at least list our cast of main characters: 



we will use the following notation: 



X 



However, 

Definition 2. ||Roj97a|| Given any w El] 

The algebraic torus (K*)^ 
The face of P with inner normal w 

The closure of the cone generated by the inner normals of 
The affine chart of Tp corresponding to the cone of Fan(P) 
The dimlP"^) - dimensional subspace o/M" parallel to P^ 
The point in Uw corresponding to the semigroup homomorphism fl Z"^' 



T 

pw 



Sp{Q) 

Div(/) 
'Dp{f,Q) 
Vp{F,V) 



where 6^ denotes the Kronecker delta 



{0, 1} mapping p 

The T -orbit of = The T -orbit corresponding to RellntP'" 
The T -invariant Weil divisor of Tp corresponding to a polytope Q with which P 
is compatible 

The Weil divisor ofTp defined by a rational function f on [K*)^ 
Div(/) + Sp{Q) = The toric effective divisor ofTp corresponding to {f,Q) 
The (nonnegative) cycle in the Chow ring ofTp defined by C(i^i1^p{fi, Pi) 

We will say that a polytope P (resp. an n-tuple V) is compatible with Q iff every 
cone of Fan(Q) is a union of cones of Fan(P) (resp. a union of cones of Fan(Pj) for each 
i) ||Kho77| , [Roj97a|| . We will also make frequent use of the natural correspondence 
between the face interiors {RellntP'"} and the T-orbits {O^} ||KSZ92| , |F\n93i |GKZ94 . 



The following lemma gives a more explicit algebraic analogy between the vertices of 
P and the maximal affine charts of Tp. 

Lemma 1. |P^oj97a|| Suppose F is a k x n polynomial system over K with support 
contained in a k-tuple of integral polytopes V := (Pi, . . . , Pk) in M". Assume further 
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that P is a rational polytope in compatible with V . Then the defining ideal in 
K[x^ leealnZ^] of U^nVp{F, V) IS {x'^fi, ... , x^'^h), for any h, . . . ,bk e 
such that h + P^"" C for allie[l..k]. ■ 

Theorems 1 and 2 will follow easily from the Vanishing Theorem for Resultants 
via following two lemmata below. 

Lemma 2. Following the notation of theorem 1, define 

VE:={Conv{E,),... ,Conv(^„)). 

Then the map x ^ x"' extends to a proper morphism (pa '■ Tpe\L{E, a) — > P]^. Also, 
ifVpj^{F,VE) has zero- dimensional intersection with (K*)^ and avoids the ambiguity 
locus, then TTa is, up to a nonzero constant multiple, the Chow form of the subvariety 
UVp,{F,Ve)) ofW],. ■ 

Lemma 3. Following the notation of theorem 1, pick points a(l), . . . , a(n — 1) G Z" 
which together generate a hyperplane in M" which is not parallel to any facet of Pe. 
Then F has infinitely many roots in {K*)^ =^ the polynomial p{u^, u_) := YYiZi '^a(i) 
is identically 0. ■ 

Sufficiently armed, let us now prove theorem 2. 

Proof of Theorem 2: Suppose that F has no roots within L{E, a) and only finitely 
many roots within (X*)". Then it follows from lemma |] that (resp. m_) divides 
TTa iff F has a root in some torus orbit with w-a>0 (resp. w ■ a < 0). More to 
the point, if cG K* then it follows similarly that m+ + CU- divides TTa iff some root 
of F lies in the closure of the hypersurface {x°' = c} in Tp^. Furthermore, lemma 1 
immediately implies that any two (distinct) hypersurfaces of this form must intersect 
precisely on the ambiguity locus. Thus tTq can not be identically and we obtain 
that M is finite. 

Now assume that F has a root in the ambiguity locus. Then, by the Vanishing 
Theorem for Resultants and what we've just learned about the hypersurface {x" = c}, 
TTa must be divisible by infinitely many distinct binomials of the form + cit_. So 
the existence of a root of F within the ambiguity locus implies that A/' = oo. Similarly, 
the existence of infinitely many roots of F within (K*)^ implies that A/'=oo as well. 
Thus we are done with the first portion of the theorem. 

As for the codimension of the space of "degenerate" F being > 2, first consider the 
F which have a root on the ambiguity locus. These F actually form a codimension 2 
subvariety. To see this, first note that F has a root in <^=^ the initial term system 
in^(F) has a root in (K*)" |[Roj974 corollary 2]. By the definition of L{E, a), we can 



then conclude via theorem 1.3 of | Stu94|| . (Although the results of ||Stu94|| are stated 



over the complex numbers, this is only a minor technicality: in this case, the results 
of ||Roj97a|| immediately imply that we can apply theorem 1.3 over any algebraically 



closed field.) 
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As for the F with infinitely many roots in {K*)^, assume temporarily that the 
coefficients of F are algebraically independent indeterminates. That the F with 
infinitely many roots in {K*)^ define a subvariety of codimension >2 must certainly 
be an old result. However, for completeness, we will give a quick proof: Following the 
notation of lemma 0, Hilbert's Irreducibility Theorem ||Lan83|| readily implies that we 



can choose a, /3, 7, 6^K* such that the polynomials p{a, (3) and ^(7, 5) are relatively 
prime in K[Ce\- So, by lemma ^ these F are indeed contained in a codimension 2 
subvariety: {Ce \ /3) =^(7, 5) =0}. ■ 

So in summary, computing tTq amounts to laminating with hypersurfaces of 
the form {c+x" = c_}. These laminae all intersect at the ambiguity locus, and toric 
infinity is precisely the union of the two degenerate laminae corresponding to (cf. 
remark ^ the ± halves of OPe- 

Proof of Theorem 1 and Corollary 1: We need only state the preceding para- 
graph more algebraically. In particular, lemma |^ has already done this for us. We 
thus obtain that when M < 00, the multiplicity of a factor c+z/+ + c_m_ of TTa is 
precisely the sum of the intersection numbers of the components of Pp^(F, Ve) lying 
in the fiber (^~^([c_ : c+]). So A/" (resp. J\f') is indeed the exact number of roots, 
counting multiplicities (resp. not counting multiplicities), of F in (K*)^. ■ 

Note that it is actually possible for a positive-dimensional component of T'p^(F, Ve) 
lying at toric infinity to make a positive (finite) contribution to A/" (or A/"'). This can 
be interpreted as an alternative way of associating an intersection number to an 
excess component | |b'ul84|| . 



6. Toric Generalized Characteristic Polynomials 

We will first comment briefiy on the combinatorial preprocessing step. The follow- 
ing example illustrates a simple fundamental case. 

Example 3 (The "Dense" Case). Suppose V is the n-tuple (rfiA, . . . ,dnA) where 
A C is the standard n-simplex and di&N for all i. It is then easily verified that 
the n-tuple D := {{O , didi} ,.. . ,{0,(i„e„}) is an irreducible fill of V ||Roj94| , ^ch94|| . 



Letting A be the vertices of A, theorem 3 then implies that 7i is a variant (over 
a general algebraically closed field) of the original GCP applied to an nxn system 
of equations with degrees di, . . . ,dn [Pan9CI|| . In particular, our F — sF* has 2n 
s-monomials, compared to n s-monomials in Canny's (/i — sxf^ , . . . , fn — sx^"). 
Note also that Conv(74) and P are homothetic and Tp = P^. Neglecting the extra 
s-monomials, setting di = l for all i, and suitably specializing the coefficients of g, we 
can then recover the usual characteristic polynomial of a matrix. 

We point out that the computational complexity of finding an irreducible fill is an 
open question. However, the connection between fills and polynomial system solving 
(not to mention specialized resultants) appears to be new and, we hope, provides 
added incentive to investigate filling. Also, even if finding a fill is difficult, this step 
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need only be done once for a given family of problems, provided E remains fixed. The 
situation where the monomial term structure of a polynomial system is fixed once and 
for all (and the coefficients may vary thousands of times) actually occurs frequently 
in many practical contexts, such as robot control or computational geometry. 

In any event, the toric GCP is an algebraic perturbation method and irreducible 
fills provide a combinatorial means of inducing "general position" into the roots of 
!Fa- For example, the following lemma implies that F* is sufficiently generic in a 
useful sense. 

Lemma 4. Following the notation and assumptions of theorem 3, for any point v 
lying in any Di, there exists a w E M"\{0} such that {i} is the unique essential 
subset of and Df = {v}. In particular, F* has exactly M.{V) roots (counting 
multiplicities) in (K*)"' . ■ 



This lemma follows easily from the techniques of [[Roj94|| , particularly section 2.5. 



Given that F* thus has maximally many isolated roots (counting multiplicities) 
and no excess components in (K*)^, theorem 3 then follows easily from the Vanishing 
Theorem for Resultants and an algebraic homotopy argument. A similar homotopy 
technique appears in | [RW96| , prloj97a| ], and the dense case (with K = C) is covered 



in | |(Jan9(J| , [Shu9^i ] . (However, toric variety language was not used in the last two 
works.) Thus, the usual GCP is (almost) the special case of the toric GCP where Tp 
is complex projective space. 

In closing, an important difference to note is that our present toric GCP is primarily 
suited for {K*)^, while the original GCP is mainly suited (in a non-sparse way) for 
affine space. To completely generalize and improve the GCP in affine space, it is 



necessary to use the affine sparse resultant and this is pursued further in |[Roj97e |. 
For instance, by replacing the sparse resultant with the affine sparse resultant, and 
using K^-counting [[Roj97a|| instead of filling, we can actually recover Canny's GCP 



in the dense case. 
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